Theory of the field-effect mobihty in amorphous organic 

transistors 



M. C. J. M. Vissenbergi'2'0 and M. Matters^'0 
^Philips Research Laboratories, 5656 A A Eindhoven, The Netherlands 
Instituut-Lorentz, University of Leiden, 2300 RA Leiden, The Netherlands 

(February 1, 2008) 

Abstract 

The field-effect mobility in an organic thin-film transistor is studied theo- 
retically. Prom a percolation model of hopping between localized states and 
a transistor model an analytic expression for the field-effect mobility is ob- 
tained. The theory is applied to describe the experiments by Brown et al. 
[Synth. Met. 88, 37 (1997)] on solution-processed amorphous organic tran- 
sistors, made from a polymer (polythienylene vinylene) and from a small 
molecule (pentacene). Good agreement is obtained, both with respect to the 

gate-voltage and the temperature dependence of the mobility. 
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Over the last decade the use of organic p-type semiconductors in field-effect transistors 
has gained considerable interest due to their potential application in low-cost integrated 
circuits. Most effort has been put into increasing the hole mobility of the semiconductor 
and increasing the on-off ratio of the field-effect transistor by optimising existing materials 
and by applying new materials. Mobilities as high as 0.7 cm^V^^s^^ and on-off ratios of 10^ 
have recently been reported in thin-film transistors of evaporated pentacene.0 Furthermore, 
attention has been focused on the improvement of the processability of these materials by 
using directly solublJ or precursor organic semiconductors.! Besides the technical applica- 
bility of organic semiconductors, their electronic and structural properties have been the 
subject of investigation as well. Interesting questions like the connection between molecular 
order and hole mobility in conjugated oligomers and polymers have been addressed. 00~i 

Experiments have indicated that the field-effect mobility of holes in organic transistors 
depends on the temperature as well as on the applied gate bias.lii This has been described by 
Horowitz et al. using a multiple trapping and release model.0 In this model the assumption 
is made that most of the charge carriers are trapped in localized states. Then the amount 
of (temporarily) released charge carriers to an extended-state transport level (the valence 
band for classical p-tjpe semiconductors) depends on the energy level of the localized states, 
the temperature, and the gate voltage. However, while extended-state transport may occur 
in highly ordered vacuum-evaporated molecular films,0 we do not expect it to play a role in 
amorphous organic films,! where the charge carriers are strongly localized. 

In the present paper, we derive a theory for the field-effect mobility in amorphous organic 
transistors, where the charge transport is governed by hopping, i.e. the thermally activated 
tunneling of carriers between localized states, rather than by the activation of carriers to a 
transport level. We use the concept of variable range hopping (VRH), i.e. a carrier may either 
hop over a small distance with a high activation energy or hop over a long distance with a low 
activation energy. The temperature dependence of the carrier transport in such a system is 
strongly dependent on the density of localized states. In a field-effect transistor, an applied 
gate voltage gives rise to the accumulation of charge in the region of the semiconducting 



layer that is close to the insulator. As these accumulated charge carriers fill the lower- lying 
states of the organic semiconductor, any additional charges in the accumulation layer will 
occupy states at relatively high energies. Consequently, these additional charges will — on 
average — require less activation energy to hop away to a neighboring site. This results in a 
higher mobility with increasing gate voltage. 

This paper is organized as follows. First, we study the infiuence of temperature and the 
infiuence of the filling of states on the conductivity of a VRH system with an exponential 
distribution of localized-state energies. Using percolation theory, we find an analytic ex- 
pression for the conductivity. This expression is then used to derive the field-effect mobility 
of the carriers when the material is applied in a transistor. Finally, our result is used to 
interpret the experimentally observed temperature and gate- voltage dependence of the field- 
effect mobility in both a pentacene and a polythienylene vinylene (PTV) organic thin-film 
transistor.! 

Let us first derive an expression for the conductivity as a function of temperature T 
and charge carrier density. At low carrier densities and low T, the transport properties are 
determined by the tail of the density of (localized) states (DOS). Our model is based on the 
following exponential DOS, 



where Nt is the number of states per unit volume, fee is Boltzmann's constant, and Tq is a 
parameter that indicates the width of the exponential distribution. We take g{e) = for 
positive values of e. We do not expect the results to be qualitatively different for a different 
choice of g{e), as long as g{e) increases strongly with e. 

Let the system be filled with charge carriers, such that a fraction 6 G [0, 1] of the 
localized states is occupied by a carrier, i.e. such that the density of charge carriers is 6Ni. 
In equilibrium, the energy distribution of the carriers is given by the Fermi-Dirac distribution 
/(e, ep), where ep is the Fermi energy (or chemical potential). For a given carrier occupation 
6, the position of the Fermi energy ep is fixed by the condition 




(1) 
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^ exp (^)r(l - T/nni + T/To), (2) 

where T{z) = dy exp {—y)y^''^ . In Eq. we have used the assumption that most 
carriers occupy the sites with energies e <^ 0, i.e., — ep ^ k-sTo. This condition is fulfilled 
when S and T are low. Note that at T = the gamma functions are unity and the carrier 
density is given by the density of states with energies lower than ep. Our approximate 
expression breaks down at temperatures T ^ Tq, where r(l — T/Tq) diverges. At such 
temperatures our assumption that only the tail of the DOS is important no longer holds, as 
the majority of the carriers is located close to e = 0. 

The transport of carriers is governed by the hopping of carriers between localized states, 
which is strongly dependent on the hopping distances as well as the energy distribution of 
the states. At low bias, the system can be described as a resistor network.i'i In this case, 
one can assign a conductance Gij = Gqc"^^^ between site i and site j where 

s^j = 2aR,, + (3) 

Here, the first term on the right-hand-side describes the tunneling process, which depends on 
the overlap of electronic wavefunctions of the sites i and j. In a lowest-order approximation, 
this tunneling process may be characterized by the distance Rij between the sites and an 
effective overlap parameter a. The second term in Eq. (|^) takes into account the activation 
energy for a hop upwards in energy and the occupational probabilities of the sites i and j. 
According to percolation theory,il0 the conductivity of the system is given by 

a = aoe-'% (4) 

where (Tq is an (unknown) prefactor and Sc is the exponent of the critical percolation con- 
ductance Gc = GqC"'^''. This Gc is determined as follows. Take a reference conductance G 
and remove all conductances with Gij < G. For high G, the remaining conductances form 
isolated clusters. These clusters increase in size with decreasing G. The critical percolation 



conductance Gc is defined as the value of G at wliicli tlie first infinite cluster is formed. This 
Gc determines the conductivity (|^), since it is the most difficult step required for transport 
through a macroscopic system. The onset of percolation, i.e., the occurrence of the first 
infinite cluster, is usually determined by calculating the average number of bonds (conduc- 
tances with Gij > G) per site in the largest cluster. This number of bonds B increases 
with decreasing G, until, at the onset of percolation, a critical number Be is reached. For a 
three-dimensional amorphous system, this number is given by B, ^ 2.8.00 The percolation 
criterion B{G = Gc) = Be can be written as 

where A^b is the density of bonds and A'^s is the density of sites in the percolating system. 
The density of bonds is given by 

Nh = J dRijdeidejg{ei)g{ej)6 (s^ - Sij) . (6) 

The density of sites Ns can be estimated by excluding all sites that cannot possibly belong 
to an infinite clusterS 

N, = J deg{e)e {scksT - |e - epl) . (7) 

Substituting Eqs. (|l|) and (||) into Eqs. (||), (^, and (0), we obtain the percolation criterion 
for our system, 

where we have assumed that the site positions are random, that most of the hopping 
takes place between tail states (— cf ^ ksTo), and that the maximum energy hop is large 
[scksT ^ ksTo). We note that our percolation criterion (^) is, up to a numerical factor, in 
agreement with previous resultsS'S where somewhat different approaches have been used 
to describe VRH in an exponential band tail. 

Combining Eqs. (0), and (^, the expression for the conductivity as a function of the 
occupation 5 and the temperature T is obtained. 



Note that the conductivity has an Arrhenius-hke temperature dependence a ~ 
exp [— i^a/l^B^)]? with an activation energy E^^ that is weakly (logarithmically) tempera- 
ture dependent. This is in strong contrast with the well-known Mott formula for VRH in 
a constant DOS, where a ~ exp — (Ti/T)^^^ .111 The temperature dependence of the Mott 
formula is a consequence of hopping over far distances and hopping to high energies being 
equally important. In an exponential DOS, however, the characteristic hop is an activated 
jump, since there are much more available states at higher energies. In fact, it has been 
shown that hopping in an exponential DOS can be effectively described in terms of activa- 
tion from the Fermi energy to a specific transport energy.liilll This explains the Arrhenius 
behavior of Eq. (H). 

In our expression (|^), the conductivity increases superlinearly with the density of carriers 
(cr ~ 5"^"/^). This is due to the filling of localized states: an increase in the carrier density 
gives rise to an increase in the average energy, thus facilitating an activated jump to the 
transport energy mentioned above. When the filling of states is not taken into account, 
i.e. when Boltzmann statistics instead of Fermi-Dirac statistics is used in Eq. (H), the 
conductivity is simply proportional to the density of carriers, m 

We now apply the obtained conductivity @ to describe the field-effect mobility /ipE 
in a transistor. In bulk material, the mobility fi of the charge carriers is given by /i = 
a{6,T)/{e6Nt). In a transistor, however, the charge density is not uniform, but it decreases 
with the distance x from the semiconductor-insulator interface. According to Eq. (^, the 
occupation 6{x) depends on the distance x through the gate-induced potential V{x), 



5(^) = 5oexp[^), (10) 



where Sq is the carrier occupation far from the semiconductor-insulator interface, where 
V{x) = 0. The variations of V{x) and 6{x) with the distance x are determined by the 
Poisson equation. For the accumulation layer, where S{x) ^ 6o, the following relation 
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between the electric field F{x) = —dV{x)/dx and S{x) can be obtained,! 

F\x) = 2k^T^N,5{x)/e,, (11) 

where eg is the dielectric constant of the semiconductor. The field -F(O) at the interface can 
be expressed in terms of the gate voltage Vg and and the insulator capacitance per unit area 
Ci through Gauss' law, 

F(0) = a^G/es. (12) 

Substituting the distance-dependent charge occupation 5{x) into Eq. (^) for the conduc- 
tivity, the source-drain current of the transistor in the linear regime {—Vd < —Vq) reads 

/ = -4^ / dxa[S (x) , T] . (13) 
L Jo 

Here, Vd is the drain voltage (the source is the ground electrode) and L, W and t are 
the length, width, and thickness of the channel, respectively. The field-effect mobility then 
follows from the transconductance (see, e.g., Ref. ^) 

_ L dl 



From Eqs. (|9|)-(p!^ we obtain the following expression for the field-effect mobility, 

...^-(^ (15) 

e \{2afB,T{l-T/To)T{l + T/To)J L2A;bTo e^J ^ ^ 

where we have assumed that the thickness t of the semiconductor layer is sufficiently large 
such that V{t) = 0. Then the field-effect mobility is independent of the thickness t as well as 
the bulk carrier occupation go- We note that the A^t dependence of the charge distribution 
in the accumulation layer is exactly cancelled by the dependence of (r{q,T). 

Let us now apply our result (|I5|) to the experimental data of Ref. 0, where the drain cur- 
rent I versus gate voltage Vg characteristics have been measured of both a pentacene and a 
polythienylene vinylene (PTV) organic thin-film transistor at a range of temperatures. The 
precursors of both organic semiconductors are spin-coated from solution on a substrate con- 
sisting of a heavily n-doped silicon (common) gate electrode, a 200 nm thick Si02 insulating 



layer (Ci = 17 nFcm~^) and a patterned gold layer as the source and drain electrodes. The 
precursors are converted into the organic semiconductors using a process described in Ref. |^. 
Typical channel widths and lengths were W = 10-20 mm and L = 2-20 fim respectively. 
The film thickness t varied from 30 to 50 nm. For both semiconductors, we use a relative 
dielectric constant = 3, which is appropriate for most organic solids. In Fig. [I| the field- 
effect mobility in a pentacene and in a PTV thin-film transistor is plotted against the inverse 
temperature for different gate voltages. Experimentally, the field-effect mobilities are deter- 
mined from Eq. ( |14|) at Vd = —2 V. The theoretical curves (solid lines) follow from Eq. (|15|), 
where we have used ctq, a, and Tq as fitting parameters. The agreement with experiment 
is quite good (the parameter values are given in Table |ID. The temperature dependence 
of /ipE, as shown in Fig. |l|, follows a simple Arrhenius behavior /xpE ~ exp [— _E'a/(^B^)], 
where the activation energy E^. depends on Vq as plotted in Fig. ||. The decrease of i?a 
with increasing (negative) gate voltage is the direct result of accumulated charges filling the 
lower-lying states. As a result, any additional charge carriers in the system will occupy sites 
with — on average — a higher energy and less energy will be required for the activated jumps 
to neighboring sites. 

The field-effect mobility in PTV is more than two orders of magnitude lower than the 
field-effect mobility in pentacene. Furthermore, the activation energy for PTV is about twice 
the activation energy for pentacene. Surprisingly, these differences cannot be attributed to 
differences in the prefactor ctq nor to the width of the energy distribution Tq, as these 
parameters have similar values for PTV as well as pentacene (see Table |ID. The main 
difference between pentacene and PTV appears to be in the overlap parameter a, which 
determines the tunneling process between different sites. We note that this key parameter is 
absent in a multiple-trapping model, where the transport is governed by thermal activation 
from traps to a conduction band and subsequent retrapping, without involving a tunneling 
step. As the length scale is smaller than the size of a molecule, one must be cautious 
not to interpret simply as the decay length of the electronic wave function. The size 
and shape of the molecules and the morphology of the organic film are expected to have 
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an important influence on the tunneling probability as well. The observed difference in 
may be due to the fact that there is more steric hindrance in the polymer PTV than in a 
system of small pentacene molecules. The better stacking properties of pentacene give rise 
to a larger area of overlap of the electronic wave functions, which results in a larger effective 
overlap in our model. 

In conclusion, we have derived an analytic expression for the fleld-effect mobility in a 
thin-fllm transistor of an amorphous organic semiconductor, using percolation theory and the 
concept of hopping in an exponential density of localized states. The calculated temperature 
dependence and gate- voltage dependence agree well with those of the the observed fleld-effect 
mobility in both a pentacene and a PTV thin-fllm transistor. According to the theory, the 
differences in the magnitude and in the temperature dependence of the fleld-effect mobility 
of pentacene and PTV transistors are mainly due to differences in the structural order of 
the organic fllms. 

We would like to thank A. R. Brown, C. P. Jarrett, and D. M. de Leeuw for providing the 
data of Fig |I|, and M. J. M. de Jong for useful discussions. Financial support from the EC 
under Esprit project 24793 Frequent and from the Dutch Science Foundation NWO/FOM 
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FIGURES 




lOOO/T (1000/K) 

FIG. 1. Field-effect mobility ^pE in a pentacene and a polythienylene vinylene (PTV) thin- film 
transistor as a function of the temperature T for different gate voltages Vg = -20 V (triangles), 
-10 V (circles) and -5 V (squares). The experimental data (symbols) are taken from Ref. |3|. The 
solid lines are according to Eq. ([l5|) , using the parameters given in Table |. 
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FIG. 2. Activation energy Eg, for the field-effect mobility in a pentacene and a polythienylene 
vinylene (PTV) thin- film transistor as a function of the gate voltage Vq. The experimental data 
(squares) are taken from Ref. ^. The solid lines are calculated from Eq. (|T5|), using the parameters 
given in Table |. 



12 



TABLES 

TABLE L The pre-exponential factor to the conductivity do, the overlap parameter and 
the width of the exponential distribution of localized states Tq for both pentacene and polythieny- 
lene vinylene (PTV) as obtained from the fit of Eq. ([l^) to the experimental data of Ref. see 
Fig-i 





ao (lO^o S/m) 


(A) 


^o(K) 


pentacene 


1.6 


2.2 


385 


PTV 


0.7 


0.8 


380 
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